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Pairing properties of even-even superheavy 7V=184 isotones are studied within the 
Skyrme-Hartree-Fock+BCS approach. In the particle-hole channel we take the Skyrme 
energy density functional SLy4, while in the particle-particle channel we employ the 
^-■^ ' seniority pairing force and zero-range (5-interactions with different forms of density de- 

pendence. We conclude that the calculated static fission trajectories weakly depend on 
the specific form of the 5-pairing interaction. We also investigate the impact of triaxiality 
on the inner fission barrier and find a rather strong Z dependence of the effect. 

o 

^ ' 1. Introduction 

In the almost fifty years since the phenomenon of superconductivity was brought 
' into nuclear structure, "'^ ^ ^ our knowledge of the nature of pairing correlati ons in 

^ . nuclei is still unsatisfactory and many fundamental questions remain (see Refs.^EE] 

and papers quoted therein). For instance, some components of the pairing inter- 
action are believed to be induced, like in solid-state superconductors!^ some are 
directly rooted in the nucleon-nucleon force. In practical calculations based on ef- 
fective interactions, one considers various pairing parametrizations. The resulting 
pairing fields have a strong influence on most low-energy properties of the nuclei^ 
and the nuclear large-amplitude collective motion. Indeed, by increasing configu- 
ration mixing and reducing the magnitude of symmetry-breaking effects, pairing 
tends to make the nuclear collective motion more adiabatic,l21 Therefore, when aim- 
ing at a quantitative understanding of fission properties of heavy and superheavy 
elements, it is important to have the pairing channel under control. 

The purpose of this study is to compare different pairing schemes that are 
currently used in the Skyrmc-Hartrce-Fock+BCS (SHF+BCS) model to describe 
nuclear superfluidity. Within this framework, we compare the seniority paring force 
and state-dependent 5- interaction (referred to as SHF-I-BCS(G) and SHF-f BCS(5), 
respectively). As in our previous paper ,^1^ we focus here on the fission properties 
of even-even isotones with 7V=184; namely, we discuss the total binding energies. 
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mass hexadecapole moments, and pairing gaps calculated along the static fission 
paths. 

2. Theoretical framework and results 
2.1. Description of pairing correlations 

In our SHF+BCS approach, we use the Skyrme energy density functional in its 
SLy4 parameterizatioJUl in the particle-hole channel, whereas two different pair- 
ing schemes were implemented in the particle-particle channel. The SHF+BCS (G) 
sch eme employs the seniority pairing force with strength parameters defined as in 
Ref^, i.e., 

G" = [19.3 - 0.084 (TV - Z)] /A , 

GP = [13.3 + 0.217 {N - Z)] /A, ^ ^ 

additionally scaled by 

G"/^' - /„/pG"/P. (2) 

In the SIIF+BCS((5) scheme, we apply the state-dependent (5-interaction '-^ with 
commonly used parameterization variants, which are summarized as 



Vg " {ri,r2) = Vq 



1 — 7] 



PO 



S{n~r2), (3) 



where po = 0.16 fm and 

{0 , for the delta interaction (DI) , volume pairing; 
1 , for the density-dependent delta interaction (DDDI), surface pairing; (4) 
i , for the mixed volume and surface pairing (MIX). 

The scaling factors of Eq. fn = 1.41 and fp = 1.13, and pairing strengths 
Vo" = 282.0 MeV, = 285.0 MeV (DI), = 842.0 MeV, = 1020.0 MeV 
(DDDI), and = 425.5 MeV, = 448.5 MeV (MIX) were adjusted to reproduce 
the experimentailS neutron (A„ = 0.696 MeV) and proton (Ap = 0.803 MeV) 
pairing gaps in ^^^Fm. As we deal with contact interactions, we use a finite pairing- 
active space defined by including fl^^P = (N or Z) lowest single-particle states for 
neutrons and protons, respectively. In the SHF-|-BCS((5) approach, the pairing gap 
is state dependent. Therefore, the average (spectral) gaps, 

were used as measures of experimental pairing gaps deduced from the odd-even 
mass staggering. In Eq. ^ Vk and Uk are the BCS occupation amplitudes (see, 

e.g., Ref.'l^' for a more detailed discussion). 

The calculations were carried out using the code hfodd 

(v.2.i9iFiwgi 

that 

solves self-consistent HF equations by using a Cartesian 3D deformed harmonic- 
oscillator finite basis. In the calculations, we took the lowest 1140 single-particle 
states for the basis. This corresponds to 17 oscillator shells at the spherical limit. 
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Fig. 1. (A) The total binding energies E^°^ (left-hand side scale) and mass hexadecapole moments 
Q40 (right-hand side scale) along the fission paths of ^^*Rfi84 calculated with the SLy4 interaction 
and four different pairing interactions: MIX, DDDI, and DI 5-interaction, and seniority pairing 
(G). (B) The neutron A„ and proton Ap pairing gaps along the fission paths shown above. 




2.2. Comparison of pairing models 

Figures [U [21 and [3] display total binding energies mass hexadecapole mo- 

ments ((34o), and neutron/proton pairing gaps (A„/p) calculated along the static 
fission paths of ^^^Rfi84, ^^^114i84, and ■^^°126i84. The fission paths were computed 
with a quadratic constrainlP^ on the mass quadrupole moment {Q2o)- Our study 
covers the oblate/prolate deformations of Q20 = —30 3 00 b wi th a step of 10 b. 

As in previous studies based on the SLy4 functionalp^^ we found that all 
nuclei considered in this work have spherical shapes in their ground states. This 
is because A'^=184 appears to be the magic neutron number in most theoretical 
models based on the Skyrme approach. ■ Due to the magic character of N^ISA 
isotones, ground-state neutron pairing gaps calculated in BCS vanish. The ground- 
state proton pairing gaps, on the other hand, show considerable variations with Z. 
They are large in the open-shell rutherfordium (Z=104), but in Z—llA the values 
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Fig. 2. Same as in Fig.[T]but now for ^^*114i84. 



of Ap are considerably reduced. As expecte d, proton pairing is much weakened for 
Z=126 which is predicted to be semi-magic.^2ll22!23' 

The calculated static fission paths have reflection-symmetric shapes, i.e., Qxo = 
for all odd multipolarities A. Furthermore, one can see that irrespective of the 
pairing interaction used, the hexadecapole moments are almost identical, and they 
gradually increase (from up to about 80 b^) along the calculated static fission 
paths. In contrast to Q40, the collective potentials (i.e., total energies E^°^ as func- 
tions of Q20) differ depending on the pairing model employed. This is particu- 
larly evident for ^*^Rfi84 (Fig. [Tl\) where the fission barrier calculated within the 
SIIF-I-BCS(G) model is significantly higher as compared to those obtained in the 
SIIF+BCS((5) variants. In the case of ^^^114i84 (Fig. OA,), the fission barriers calcu- 
lated with the SIIF-I-BCS(G) and SIIF+BCS((5) interactions are considerably closer 
to one another than in the case of ^^*Rfi84. Furthermore, for ■^^*'l26i84 (Fig.[5J'V), all 
fission barriers calculated with both pairing models are almost identical. This can 
be attributed to the behavior of proton pairing along the fission paths. Indeed, in 
^®*Rfi84 there is a large systematic difference between Ap values in SIIF-I-BCS(G) 
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Fig. 3. Same as in Fig.[T]but now for 3iOl26i84. 

and SHF+BCS((5) variants, with the seniority-pairing model producing considerably 
larger pairing gaps. This difference decreases when going towards ^^"126184 in which 
An/p obtained within the SHF-I-BCS(G) model are much closer to those obtained 
within the SHF+BCS((5) model. This result indicates that the isospin dependence 
of seniority pairing strengths given by Eq. ([T]) is too strong and thus unrealistic. 
Another interesting observation is that neutron/proton pairing gaps (hence, cor- 
responding potential energies) calculated within the SHF-|-BCS(5) framework are 
very similar, regardless of the parameterization variant used (MIX, DDDI or DI), 
see Eq. (gl). 

2.3. Fission barriers for the even-even N —184 isotones with 
S -pairing interaction 

Figured] compares the total binding energies (_E*°*) and the mass hexadecapole mo- 
ments {Q40) calculated along the fission paths for twelve even-even A^=184 isotones. 
Here we use the MIX parameterization of the (5-interaction. We found that all of 
the super-heavy elements (SHE) studied here have reflection-symmetric static fis- 
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Fig. 4. The total binding energies E^°^ (• - scales on the left-hand side) and mass hexadecapole 
moments Q40 (♦ - scales on the right-hand side) along the fission paths, calculated with the SLy4 
interaction and the MIX parameterization of the (5-interaction for the even-even Af=184 isotones 
of 288<A<310. The differences between open and solid symbols illustrate the effects of triaxiality 
on the inner barriers. 



sion paths and are spherical in the ground states, and the Q4Q moments calculated 
along the static fission paths follow the same pattern; i.e., their values continuously 
increase from up to about 80 b^. 

All collective potentials shown in Fig. [4] have two-humped shapes and similar 
widths. Only in the case of ^*^Rfi84, can one see an additional small third external 
barrier. The outer barrier is systematically reduced with Z. The differences between 
open and solid symbols illustrate the effect of triaxiality on the inner barriers. The 
strongest effect is predicted for the nucleus ^"'^''I26i84, where the barrier is lowered 
by more than 3MeV due to triaxiality. However, for lighter isotones with Z < 114, 
the role of triaxiality is modest. 

For completeness, it is interesting to look at the neutron (A„) and proton (Ap) 
spectral pairing gaps calculated along the static fission paths shown in Fig.[31 They 
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Fig. 5. The total binding energies (•j scales on the left-hand side) and the spectral neutron 

(A„) (a) and proton (Ap) (▼) pairing gaps (scales on the right-hand side) for the Af=184 isotones 
shown in Fig. U] 



are shown in Fig. [5l As already mentioned, the ground-state neutron pairing gaps 
(A„) (a) are equal to in all the N—ISA isotones, whereas the proton gaps (Ap) 
(T) change from 1.3 MeV in '^^^Rhsi) to less than 0.5 MeV in 288126134. In the 
sub-barrier regions of the fission paths (shadowed regions in the plot), neutron and 
proton spectral pairing gaps fiuctuate around 1 MeV. 

The size of a fission barrier is a measure of stability of the nucleus against spon- 
taneous fission. The static fission barriers shown in Fig. [1] suggest a possible increase 
of the stability against spontaneous fission for 118 < Z < 124. This predicted in- 
crease of the stability supports the results obtained in our previous paper ,^1^ where 
the same nuclei were considered within the SIIF+BCS(G) framework, but with the 
pairing strengths of Eq. |^ scaled to reproduce the pairing gaps of the finite-range 
droplet model (FRDM).l2il 
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3. Conclusion 

The purpose of this work was to study the sensitivity of static fission trajectories 
in A''=184 isotones to the choice of pairing interaction in the SHF+BCS model. We 
found that the results (energies, hexadecapole moments, pairing gaps) arc fairly 
insensitive to the assumed density dependence of the pairing (5-interaction (DI, 
DDDI, or MIX). On the other hand, we found rather large differences between 
results obtained in (5-pairing and seniority-pairing variants. This suggests that the 
isospin dependence of the seniority-pairing strength used is unrealistic. In general, 
it is always advantageous to consider a more realistic interaction in the pairing 
channel (such as the density-dependent 5 force) which does not depend explicitly 
on mass number and neutron excess. 

In the second part of the paper, we examined static fission paths calculated 
within the SHF-I-BCS(^) framework and MIX pairing parameterization for twelve 
even-even A''=184 isotones. We found that the reduction of fission barriers due 
to the appearance of triaxial deformation strongly depends on Z. Furthermore, 
the analysis of the static fission barriers suggests the increased stability against 
spontaneous fission for 118 < Z < 124. 
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